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ABSTRACT 

We report 60 and 90 /jm observations of 7 millisecond 
pulsars with the ISOPHOT instrument and describe our 
analysis procedures. No pulsars were detected, and typi- 
cal (3(7 ) upper limits are 150 mJy. We combine our results 
with others in the literature and use them to place con- 
straints on the existence of protoplanetary or dust disks 
around millisecond pulsars. 
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1. INTRODUCTION 
The first extrasolar planets discovered were found around 



PSR B1257+12 



the m illisecond pulsar PSR B1257+12 (Wolszczan fc Frail 



1992 ). Vari ous mechanisms have b een proposed for their 



formation ( Phinney fc Hansen 1992 ) , but all generally rely 
on an accretion disk around the pulsar within which the 
planets form. Millisecond pulsars themselves are under- 
stood to be neutron stars that have been "spun-up" by 



the a ccretion of matter from a companion ( van den Heuvel 
|1995| ) 



Both of these arguments suggest that dust disks may 
exist around millisecond pulsars. These dust disks may 
exist regardless of whether or not the pulsar is orbited by 
planets. The planetary formation process is unlikely to be 
100% efficient, and a dust disk may represent the remnant 
debris from the formation of the planetary system. Alter- 
nately, not all accretion disks may give rise to planetary 
systems, and a dusk disk would represent the remnant of 
the initial accretion that spun-up the pulsar. 

A modest number of unsuccessful searches for infrared 
emission from dust around millisecond pulsars have been 
conducted. Figure [l] compares the current observational 
limits for PSR B1257+12 with the predicted emission lev- 
els from a dust disk model. The limits for other pulsars 
are similar. Figure |l| illustrates the need for far-infrared 
observations. 

2. OBSERVATIONS 

In order to generate a millisecond pulsar ISOPHOT tar- 
get list, we compiled a list of millisecond pulsars known 
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Figure 1. Upper limits on the infrared and sub-millimeter 
emission from PSR B1257+ 12. Limits are from Foster 
& Fischer (1996, FF96 and IRAS), Phillips & Chan- 
dler (1994, PC94), and Greaves & Holland (2000, GH00). 
The solid line shows the expected emission from a 300 M§ 
dust disk composed of dust particles 0.1 [im in size and 
heated by 1% of the spin-down luminosity of the pulsar 
(Foster & Fischer 1996). The pulsar is assumed to be at 
a distance of 620 pc. 



prior to 1994 August and with dis tances less than 1 kpc 
Distances are estimated from the Taylor fc Cordcs 1993 



model and should be accurate to approximately 25%. Most 
of these millisecond pulsars lie at high Galactic latitudes. 

Of these, seven were observed with the ISOPHOT in- 



strument (Lemke et al. 1996) onboard the ISO satellite 
flKesslcr ct al. 1996Q . Table |l| summarizes the observing de- 
tails. All of the observations used the P32 observing mode 
with the C100 detector. In this mode the spacecraft was 
commanded to cover a series of raster pointings around 
the nominal pulsar position. At each raster pointing an 
internal chopper pointed the beam toward 13 adjacent sky 
positions. The throw of the chopper was larger than the 
offset between raster pointings. The result was that, in 
general, an individual sky position within the raster was 
observed multiple times or oversampled. Before and after 
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each observation of a pulsar, an internal calibration source 
was observed. 



Table 1. Pulsars Observed 



Name 
PSR 



A 



D 

(PC) 



(Le) 



P32 
Raster Time 

00 



J0034-0534 


90 


1000 


10 


3 x 


8 


1402 


J1640+2224 


60 


1190 


0.88 


3 x 


6 


1012 




90 


1190 


0.88 


3 x 


8 


848 


J1730-2304 


60 


510 


< 0.35 


3 x 


6 


1590 




90 


510 


< 0.35 


3 x 


8 


1232 


B1855+09 


60 


900 


1.1 


3 x 


6 


1012 


J2124-3358 


60 


240 


1.7 


3 x 


6 


1012 




90 


240 


1.7 


3 x 


8 


848 


J2145-0750 


60 


500 


< 0.048 


3 x 


6 


1590 




90 


500 


< 0.048 


3 x 


8 


848 


J2322+2057 


60 


780 


0.62 


3 x 


6 


1804 




90 


780 


0.62 


3 x 


8 


1402 



The analysis of the pulsar observations largely followed 
the standard ISOPHOT analysis. The key difference was 
the amount of "deglitching" performed. Glitches result 
from cosmic rays striking the detector or secondary elec- 
trons produced by spacecraft materials struck by primary 
cosmic rays. Failure to remove glitches can corrupt the re- 
sponsivity drift correction of all data, not just those con- 
taining the glitches themselves. The standard ISOPHOT 
analysis pipeline removes glitches but does so without 
making use of the redundancy implicit in the oversampled 
P32 observations. 

Deglitching proceeded in the following fashion. Within 
each spacecraft pointing the chopper would sweep past 
a particular sky position multiple times (typically 3-5 
times). For each sky position, the median signal level was 
determined, then subtracted from all observations at that 
sky position. The observations from all sky positions were 
then combined to form a signal strength histogram. A sig- 
nal strength threshold was specified, and signals above 
this level were eliminated. Typically 3%-10% of the signals 
were eliminated in this stage. Depending upon the number 
of chopper sweeps per spacecraft pointing and deglitching 
prior to this stage, the median signal strength per sky po- 
sition could not always be determined accurately. In these 
cases, additional manual deglitching was done to remove 
any remaining outlier signals. Our use of the observations 
of the internal calibration sources followed the standard 
ISOPHOT analysis pipeline. 

After deglitching and calibration using the internal cal- 
ibration sources, mapping was done within the ISOPHOT 
Interactive Analysis package. Measurements from the indi- 
vidual detector pixels were co-added to form a sky image, 



with the contributions from the individual detector pixels 
weighted by their distances from the image pixels. Doing 
so takes into account the beam profile falling on each de- 
tector pixel. We also employed a median flat field, which 
has the effect of reducing substantially our sensitivity to 
any extended emission in the field but increasing our sen- 
sitivity to point sources, as desired for this program. 

As an example of the resulting images produced, Fig- 
ure U shows the image of PSR B1855+09. In no case have 
we identified a point source at the location of a pulsar. Uti- 
lizing the inner quarter of the image, we determined the 
rms noise level. We take our upper limits to be 3 times 
this rms noise level. Typical (lcr) values are 50 mJy, with 
higher values seen for pulsars close to the Galactic plane. 




RIGHT ASCENSION (J2000) 

Figure 2. PSR B1855+09 at 60 fim. The rms noise level 
is 50 mJy, and the contour levels are 1.5, 1.6, 1.7, 
and 1.8 Jy. The cross marks the location of the pulsar, 
and the size of the cross is 5000 times larger than the un- 
certainty in the pulsar's location and proper motion. 



3. RESULTS 

Foster & Fischer (1996) developed a model for the in- 
frared emission from a dust disk around a millisecond 
pulsar. Their model assumes that the disk consists of par- 
ticles of a uniform radius a heated by a fraction / s d of 
the pulsar's spin-down luminosity L s( j. The total mass of 
the disk is m^. While the model is simplistic — an actual 
dust disk presumably consists of particles with a range of 
sizes, the heating mechanism is left unspecified, and non- 
equilibrium effects such as stochastic heating are ignored — 
we believe that this simplicity is justified given the uncer- 
tainties of the heating mechanism and of the environs of 
a millisecond pulsar. 

For / s d ~ 1%, typical dust temperatures are predicted 
to be T « 10-50 K for disks having ma ~ 100 and a ~ 
1 fim and illuminated by a pulsar with L s( j ~ 1 L Q . These 
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temperatures are considerably lower than those assumed 
(T « 150 K) by Phillips & Chandler (1994), who esti- 
mated disk temperatures by scaling from observations of 
T Tauri stars. The lower temperatures result from our as- 
sumption of a weaker coupling between the pulsar's spin- 
down luminosity and the disk. Phillips & Chandler (1994) 
considered disk temperature to be a major uncertainty in 
converting from measured flux densities to inferred disk 
masses. Accordingly, our assumption of a weaker coupling 
means that larger disk masses can be tolerated without 
violating the observational constraints. 

In addition to these observations with ISO, other in- 
frared and sub-millimeter observations of millisecond pul- 
sars have been conducted. Those observations most rel- 
evant to our sample of millisecond pulsars are those by 
Foster & Fischer (1996) at 10 /im and Greaves & Hol- 
land (2000) at 850 /im. Unfortunately, there is little over- 
lap between these three samples of pulsars. Most of the 
pulsars that have been observed between 10 and 850 /im 
have been observed at only one or two wavelengths. In 
general, it is therefore not possible to constrain all three 
parameters of this model with the existing observations. 

We therefore adopt an approach in which we infer 
limits on two parameters of the Foster & Fischer (1996) 
model for fiducial values of the third parameter. Else- 
where we shall consider the detailed implications for the 
well-studied objects from our sample ( Lazio et al. 2001 ). 
Here, as an example, we consider the millisecond pulsar 
PSR J0034-0534 jBailes et al. 1994[ ) which has a probable 
white dwarf companion, is at a distance of 1 kpc, and has a 
spin-down luminosity of 10 Lq. Greaves & Holland (2000) 
placed a 2a limit of 3.7 mJy at 850 /im, and we place 
a 2cr limit of 50 mJy at 90 /im. Figure || shows the al- 
lowed region of the disk mass-grain size plane given these 
observational limits and an assumed heating efficiency of 

/ad - !%• 

Allowed regions in the md-a plane occur for one of two 
possible reasons. First, the peak of the dust disk emission 
may appear shortward of 90 /im, where no constraints 
exist for this pulsar, with the Rayleigh-Jeans tail of the 
emission falling below the two measured values. This re- 
gion is to the lower left in Figure Second, the peak of 
the emission may appear between 90 /im and 850 /im, but 
with a magnitude comparable to that measured at 90 /im 
so that the Rayleigh-Jeans tail again does not violate the 
850 /im limit while the Wien tail of the emission does not 
violate the 90 /im limit. This region is to the lower right 
in Figure ||. Obviously, a lower value of f s d would produce 
larger allowed regions in the md-a plane. 
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Figure 3. Allowed disk masses and grain sizes for a dust 
disk orbiting PSR J0034— 0534- The cross-hatched region 
indicates disk masses and grain sizes that do not violate 
the observational limits at 90 and 850 [im. A heating effi- 
ciency of / s d = 1 % has been assumed. 

lied on considerably more "deglitching" than is usually 
the case. No pulsars were detected, and typical (3c) up- 
per limits are 150 mJy. 

The current set of pulsars for which mid-infrared, far- 
infrared, and sub-millimctcr observations exists is dom- 
inated by pulsars having measurements at only 1 or 2 
wavelengths. Measurements at additional wavelengths are 
required in order to use existing dust disk emission models 
to place meaningful constraints on the presence or absence 
of a circumpulsar dust disk. 
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4. CONCLUSIONS 

We have reported 60 and 90 /im observations of 7 mil- 
lisecond pulsars with the ISOPHOT instrument. We have 
described our analysis procedures, which utilized the stan- 
dard ISOPHOT Interactive Analysis package but also re- 
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